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We construct an equation of state including the hadron-quark phase transition. The mixed phase 
is obtained by the Gibbs conditions for finite temperature. We adopt the equation of state based on 
the relativistic mean field theory for the hadronic phase taking into account pions. As for the quark 
phase, the MIT bag model of the deconfined 3-flavor strange quark matter is used. As a result, 
qq ■ our equation of state is thermodynamically stable and exhibits qualitatively the desired properties 

of hadron-quark mixed matter, such as the temperature dependence of the transition density. The 
f^} pions raise the transition density because they make the equation of state softer. Using the equation 

£SJ ■ of state constructed here, we study its astrophysical implications. The maximum mass of compact 

stars is investigated, and our equation of state is consistent with recent observations. We also 
compute the collapse of a massive star with 100 solar masses (Mg) using our equation of state and 
, find that the interval time from the bounce to the black hole formation becomes shorter for the 

model with pions and quarks. The pions and quarks affect the total energy of the emitted neutrinos 
because the duration time of the neutrino emission becomes shorter. The neutrino luminosity rises 
under the effect of pions since the density of the proto-neutron star becomes high. 
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$H ■ I. INTRODUCTION 
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, ^ , ' It has been theoretically suggested that hadronic matter undergoes a deconfinement transition to quark matter at 
high temperature and/or high density. It is also inferred that the transition occurs inside compact stars [H, @, 0). 
Stars with a quark matter central region and a hadronic matter mantle are called hybrid stars. While the equation 
of state (EOS) for the hadron-quark mixedphase is not yet fully understood, structures and maximum masses of the 
hybrid stars have been well studied 0, [f| [a 0, H, M, floL lllj- for review Ref. [HI- For the phase transition of a single 
substance, like a liquid-vapor transition of H2O, it is simple to treat the mixed phase. On the other hand, in the 
hadron-quark transition, the substance is composed not only of u quarks but also of d quarks [![• Incidentally, some 
authors use the simple Maxwell construction for the studies of hybrid stars 0, ■ In this scheme, Gibbs conditions for 
two components are not satisfied completely and EOS's of two phases are connected from the plot of pressure versus 
chemical potential. 

This transition may play an important role in the gravitational collapse of a star, such as a core collapse supernova. 
In addition, there is a possibility of forming a black hole by the transition because the compact stars have maximum 
masses. In this case, the EOS with finite temperature and including neutrinos is needed whereas the EOS can be 
calculated under zero temperature and neutrino-less /3 equilibrium for hybrid stars. The inclusion of neutrinos is 
needed because they will be trapped at a very high density as the hadron-quark transition occurs. Moreover, for the 
gravitational collapse of a star, the EOS should be constructed in a consistent method in the wide range of density. 
Recently, dynamical simulations of the stellar core collapse including the transition have been done: however, their 
EOS is obtained at zero temperature and they adopt the Maxwell construction for mixed phase [HI, [l4j]. So far, the 
stellar collapse including the hadron-quark phase transition with a finite temperature EOS and neutrinos has not 
been studied. 

A study on the mixed phase and hot hybrid stars using Gibbs conditions for finite temperature with neutrinos 
exists; however, it does not compute the dynamics of the stellar collapse and its EOS does not include low density 
regime where heavy nuclei appear || . In Ref. [f| the relativistic mean field model of the Walecka type is adopted for 
the EOS of hadronic matter. This model does not take into account the non-linear self-coupling terms of a meson and 
u> mesonjwhich are essential to reproduce the properties of nuclei quantitatively and of dense matter in a reasonable 
manner [151 ]. In addition, pions are not included in its hadronic phase while they may affect the transition density. 
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The color dielectric model is used to describe quark matter. We will revisit this model in Sec. IIII Al comparing with 
our model. 

In this study, we construct the EOS of hadron-quark mixed matter including neut rinos for fin ite temperature and 
perform the computations of the stellar collapse and black hole formation. An EOS bv lShen et al.l (Shen EOS) fl6l. []~7j] 
and the MIT bag model [l8( are adopted for nuclear matter and quark matter, respectively. Shen EOS is based on the 
relativistic mean field theory which includes the non-linear self-coupling terms and takes into account non-uniform 
matter. It is noted that the table of EOS covers a wide density and temperature range and has been often used in 
astrophysical simulations. Contributions of thermal pions are added to the original Shen EOS in the current study 
to examine its effect in the minimum model which assumes that an effective mass of pions is equal to their rest mass 
in vacuum. Incidentally, hyperons are not included in our hadronic EOS, however, we are planning to investigate 
their effects in future work 20]. The MIT bag model is characterized by a parameter B, which is called the bag 
constant and related to the degree of interaction. We calculate the EOS for several values for B and investigate its 
dependence. We show that our EOS is thermodynamically stable because Gibbs conditions for each particle species 
are completely satisfied for the mixed phase. After examining the neutron star properties, we apply our EOS to the 
numerical simulation of the stellar core collapse. It is noted that neutrinos are treated consistently in our simulation. 
Thus we can evaluate the total energy of the neutrinos emitted from the collapse. We show that the hadron-quark 
transition affects the dynamics and neutrino emission of the collapse. 

In this paper, we perform the computations only for one progenitor model and limited cases of EOS's are used for 
the hadronic phase and quark phase. This paper is the first stage of the project and we are preparing to study the 
progenitor dependence and EOS dependence [20| of the collapse and neutrino emission. The final goal of our 
study is to investigate the effects of the phase transition on the stellar collapse and the possibilities to p robe the EOS 
of hot dense matter from the emitted neutrinos, as demonstrated for pure hadronic matter in Refs. |2U [22|. 

The organization of this paper is as follows. We devote Sec. [II] to the formulation of the EOS for the hadron-quark 
mixed phase. We also introduce EOS's of the pure hadronic and quark matter and present treatments of muons. In 
Sec. IIIH we assess our EOS and discuss the dependence on the bag constant, B. As astrophysical implications, we 
evaluate the maximum mass of the compact stars and investigate the effects of hadron-quark phase transition on the 
stellar core collapse. Finally, a summary is given in Sec. IIV1 

II. FORMULATIONS 

In this section, we construct the EOS for hadron, quark and their mixture as a function of baryon mass density ps, 
electron fraction Y e and temperature T. The EOS calculated in this study covers 10 5 ' 1 g cm~ 3 < ps < 10 17 g cm -3 , 
< Y e < 0.56 and MeV < T < 400 MeV. The lower limit of ps and the range of Y e are the same as those of the 
Shen EOS, which is adopted for the pure hadronic matter in our model. We assume that the pure hadronic matter 
exists under the transition density and the pure quark matter exists for a higher density than that of the end point 
of the mixed phase. In the following, after we introduce the EOS for pure hadronic matter and pure quark matter, 
we give the formulations of the EOS for mixed matter. Finally, we also describe the treatment of muons. 

A. EOS for pure hadronic matter 

We adopt the Shen EOS f° r the pure hadronic matter. This EOS is based on the relativistic mean field theory 

and the table of EOS covers a wide baryon-mass density range, pb — 10 5 - 10 15 4 g cm" 3 , and temperature range, 
T = - 100 MeV. An inhomogeneity of the matter is also taken into account using the Thomas-Fermi approximation 
because heavy nuclei appear in T < 15 MeV and densities below the nuclear matter density. It is noted that the 
density where the hadron-quark transition occurs is higher than the nuclear matter density for T < 15 MeV in our 
model. Thus the hadronic matter at the transition point is composed of dissociated protons and neutrons. It is also 
noted that we extend the EOS to T > 100 MeV by a method fully consistent with the Shen EOS performing the 
calculations based on the relativistic mean field theory. In this temperature regime, we can regard that the nuclear 
matter is uniform and the scheme for computations is given in Refs. [23], Hi] in detail. The high temperature regime 
will be revisited in Sec. IIII Al 

We add the effects of charged pions, 7r ± , and neutral pions, ir° , to the EOS for the pure hadronic matter. The 
original Shen EOS is based on the mean field theory without the pions. It is noted that the pion field vanishes in 
the mean field approximation because it is a pseudoscalar particle [251 ]. In the realistic dense matter, pions may exist 
though large uncertainties still exist. Thus we examine their qualitative effects by treating them in the minimum 
model. Here we employ the method in Refs. H^, [2(| under the assumption that an effective mass of pions is equal to 
their rest mass in vacuum. Strictly speaking, pions feel an repulsive potential in the nucleons and their effective mass 
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gets larger than that in vacuum. In this case, the pion population is suppressed. Thus our assumption corresponds to 
one extreme case and the situation without pions does the other extreme case [25], [2(| . The realistic condition should 
be between these two cases. 

The pions are handled as an ideal boson gas and the chemical potentials of the neutral and charged pions are 
determined as follows. The chemical potential of the neutral pion is always p„o = because it is a self-conjugate 
particle and created by the pair process of two photons. As for the chemical potentials of the charged pions, we 
set them as p n - = p n — p P and p n + = — fi v - , where p p and p n are the chemical potentials of proton and neutron, 
respectively. When the charge chemical potential, p n — p, p , exceeds the rest mass of charged pions, m n ±, at high 
densities, the charged pions are condensed and the charge chemical potential decreases to their rest mass. In fact, 
there is such a regime in the original Shen EOS table. It is well known that this threshold is sensitive to the density 
dependence of the symmetry energy. We remark that the density dependence is strong in the relativistic nuclear 
many body frameworks such as the relativistic mean field theory for Shen EOS, as compared with the non-relativistic 
counter part [28| . 

In the following, we replace the original Shen EOS given as a function of baryon mass density ps, proton fraction 
Y p and temperature T with our hadronic EOS given as a function of ps, Yc and T. Yc is a charge fraction and 
defined as 

Y C = ^L^EZ (1) 
n B 

where hb and n p are the baryon number density and the number density of protons, respectively, n„- is a net number 
density of the charged pions, which is the difference of the number density of 7r~ to that of 7r + , and calculated from 
the Bosc-Einstcin distribution function: 
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Amp dp, 



{exp (y / m2 ±c 4 +p 2 c 2_ / ^^ /(fc B T)} - 1 {cxp (^±0* + pV + M7r -) /(AfeT)} 

(2) 

where h, fee and c are the Planck constant, the Boltzmann constant and the velocity of light, respectively. When the 
pions are not condensed, we find p^- and the corresponding state of the nucleons in the Shen EOS for given pb, Yc 
and T. In this process, we determine the relations of n p , p p and p n from the Shen EOS. Using fj,„- found above, the 
pressure and energy density of the pions are calculated. 

On the other hand, in case of the pion condensation, we fix the chemical potential as p^- = m^±c 2 and calculate the 
number density of the "thermal" pions, n th _ , in Eq. The pressure and energy density of "thermal" pions are given 
from the Bose-Einstein distribution function as in the case without the pion condensation. At the same time, we can 
determine the EOS of nucleons for given ps andTusing the Shen EOS under the condition of p n —^p = m 7! ±c 2 {= p, n -). 
Having Y p and n p fixed, we can get the net number density of the pions, for given Yc using Eq. ([T]). Here we 
set the number density of the "condensed" pions as n cond = n^- — n th _ . These condensed pions contribute not to 
the pressure but to the energy density by their rest masses. We can determine the EOS including condensed pions 
approximately. 

Incidentally, when the electron-type neutrinos are trapped and in equilibrium with other particles, their chemical 
potential should be given as Eq. (|8ap . which is expressed later. For high density regime where the hadron-quark 
transition occurs, neutrinos are fully trapped and we show the results including trapped neutrinos in Sec. IIII Al On 
the other hand, since neutrinos are not trapped at least for the onset of gravitational collapse, we follow the time 
evolution of neutrino distributions through neutrino reactions for numerical simulations in Sec. IIII CI In addition to 
the quark degree of freedom, hyperons and kaons also may be important in the high density regime. In particular, 
since hyperons will appear before the hadron-quark phase transition [27I . |23 |. we are investigating their effects in 
future work 12011. 



B. EOS for pure quark matter 

We adopt the MIT bag model [l8[ as the pure quark matter. This is a phenomenological model which describes 
the nature of the confinement and the asymptotic freedom of quarks. In this model, free quarks are confined in the 
"bag" and this "bag" has a positive potential energy per unit volume [18|. Thus the thermodynamical potential is 
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expressed as £Iq = Qq + BV, 1 where V and fio are the volume of the "bag" and the thermodynamical potential of 
free quarks as ideal fermions, respectively. The bag constant, B, is a parameter characterizing this model, and we 
investigate its dependence later. From thermodynamical relations, we can calculate the number density uq, pressure 
Pq and energy density eq for the quark matter with the temperature T as 

n Q = T,^J o {F?(l>)-Ff(p))4irp 2 dp, (3a) 

f h Jo sJmjc^+p^c 2 V ' 1 1 

£ Q = E |s I \j™Y+V 2 c 2 (F+ (p) + FJ (p)) 4Vdp + B, (3c) 

where Fj(p) and FJ(p) represent the Fermi-Dirac distribution functions for particle and antiparticle, respectively, 
and they are expressed as 

Ff{ P ) = - ( , , 1 , , • (4) 

[exp [^m)c A +pi C 2 TM/J AW} + 1 

The subscript / denotes the flavor of quarks, and we take into account three flavors, namely, u-, d- and s quarks. The 
statistical weight is g = 2 x 3. nif and [if are the mass and the chemical potential of / quark, respectively, and we 
adopt m u c 2 = 2.5 MeV, m d c 2 = 5 MeV and m s c 2 = 100 MeV in this study (30j . 

As mentioned already, the EOS for the pure hadronic matter is given as the function of the baryon mass density 
Pb, the charge fraction Yq and the temperature T. We can rewrite the EOS of the MIT bag model as the function of 
these three independent variables. For convenience, we define the "baryon" number density and the "baryon" mass 
density of the quark matter as 

n B,Q = — , (5a) 
Pb,q = ™ unit — > (5b) 



where the value of the atomic mass unit m un it = 931.49432 MeV is used as in the Shen EOS. Since the electric charges 

2 _ 1 

3' 3 



of u-, d- and s quarks are |, — i and — i of that of a proton, respectively, the charge fraction of the quark matter is 



as follows: 

2n„ - rid - n s 2n u ~n d -n s 

Y c,Q = = , 6 

n u +n d + n s n Q 

where n/ is a number density of / quarks. When the f3 equilibrium is satisfied, p d is equal to (j, s , which is expressed 
later as Eq. (|8d[l . Thus we can determine fi u , fi d and p s for given pb.q, Yc,q an d T, and we can express Pq 
and eq as functions of pb,q, Yc,q and T. We can also calculate other variables, such as the entropy density sq = 
Pq ^ £q + ij; J2f n ft 1 f' the Helmholtz free energy per unit volume Tq — eq — Tsq and the Gibbs free energy per unit 
volume Qq = Tq + Pq. We will use these variables later in this paper. 



C. EOS for hadron-quark mixed phase 

Following Ref. 0, we first show the conditions of the equilibrium between phases in the heat bath with the temper- 
ature T. Here, we deal with hadronic matter, quark matter, electrons, electron type neutrinos and mu-type leptons, 
whose treatment is stated later in Sec. Ill Dt We assume that the equilibrium is achieved not only by the strong 
interactions but also by the weak interactions. This is because only u- and d quarks are deconfined from protons and 



1 Some authors use a formula, which includes the lowest-order gluon interaction with the coupling constant a a [ij |13|. [i"4|. However, the 
validity is not guaranteed for the densities in which we are interested because the perturbative approach is valid only in the high energy 
limit. Hence, we do not take into account these effects in this study. 
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neutrons and s quarks are created by the weak interactions. Moreover, we assume that the neutrinos are completely 
trapped owing to high density. In other words, the diffusion time scale is much longer than the reaction time scale. 
Thus the following reactions are in chemical equilibrium: 

p + e~ < — > n + v e , (7a) 

n < — > u + 2d, (7b) 

p < — ► 2u + d, (7c) 

d < — > u + e~ + v e , (7d) 

s < — ► u + e~ + v e , (7e) 

u + d < — ► it + s, (7f) 

and the relations of the chemical potentials are given as 

Pp + Pe = Pn + Pu, (8a) 
Pn = Pu + 2 Md, (8b) 

= 2 /"« + Md. (8c) 
Pd = Ps- (8d) 

It is noted that electrons reside in both phases and the chemical potential for each phase coincides with each other. 
This is also true for neutrinos. It goes without saying that two phases are also in mechanical equilibrium. For 
simplicity, we ignore the surface tension and the screening of the charged particles though they may affect the EOS 
[3lj |. Thus we require the condition, 

Ph=Pq, (9) 

where Ph and Pq are the pressures of hadronic and quark phases, respectively. It is noted that we do not need to 
take into account the contribution of leptons in Eq. ([9]), because their chemical potentials and pressures are the same 
for each phase. 

Next, we relate the independent variables in each phase to the average values using a volume fraction of the quark 
phase, x- For instance, a pure hadronic phase corresponds to \ — and a pure quark phase to x = 1. Then the 
baryon mass density of the mixed phase, ps, is written as 

Pb = 0--x)Pb,h + XPb,q, (10) 

where pb,h is the baryon mass density of the hadronic phase and pb,q is that of the quark phase and given in Eq. 
(|5bp . We note that, the charge fractions of both phases, Yc,h and Yc,Q, do not necessarily coincide with each other. 
Since the mixed phase is charge neutral as a whole, we can relate these values with the electron fraction of the mixed 
phase, Y e , as 

Y e p B = (1 - x) y c,hPb,h + X y c,qPb,q- (11) 

It is noted that the charge neutrality is not required for each phase independently [J]. This fact gives an essential 
difference from the phase transition of a single substance like a liquid- vapor transition of H2O. For instance, the 
pressure of the hadron-quark mixed phase is not constant in an isothermal process. Now the construction of the EOS 
for hadron-quark mixed phase is reduced to the determination of the volume fraction, \, for given ps, Y e and T. 
Mathematically, it is equivalent to solving the system of 5 equations (|8bl I5cl QUI fTT|) for 5 variables, namely, x> 
Pb.h, Yc\h, Pb.q and Yc,q- Here Pb and Y e are given, p n , p p and Pjj are the functions of pb.h and Yc,Hi and p u , 
Pd and Pq are the functions of pb,q and Yq,q- 



D. Treatment of muons 



As mentioned already, the treatment of muons and their anti-particles is important because their compositions 
affect the EOS. In this study, we examine two situations. 

The first one corresponds to the EOS for stellar core collapse. Here, we call it supernova (SN) matter. In the case 
of an ordinary supernova, it is known that not only the net muon fraction but also the muon and anti-muon fractions 
themselves are minor because the rest mass of the muon is large (m M = 105.66 MeV). Then, they are omitted in 
simulations for the ordinary supernova. On the other hand, we may not be able to neglect their effects for EOS 
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in case of black hole formation, because the density and temperature would be higher than those for the ordinary 
supernova [32l. l33l|. 

During the core collapse and bounce, muons behave as follows. Since the charged current reactions for mu-type 
leptons are not efficient before neutrino trapping, mu-type neutrinos (y^) and their anti-particles (v^) are produced 
mainly by pair processes, and their reactions are almost identical. Thus one can safely posit that the net muon 
fraction is zero during the initial collapse before neutrino trapping. When neutrinos are fully trapped, the charged 
current reactions may become efficient. Therefore, muons and mu-type neutrinos are in chemical equilibrium, and 
Dfj, cannot be regarded as identical, and the net muon fraction can become non-zero. 

In this study, in order to add the effects of muon and anti-muon pairs in a tractable manner, we do not assume 
chemical equilibrium for muon type leptons 2 but suppose that the net muon fraction is zero (i.e., [i^ = 0). The 
approximation /i M = is reasonable because the pair population dominates in higher temperature regime, where we 
cannot neglect the effects of muon-type leptons for EOS in the black hole formation. Incidentally, this assumption is 
consistent with the methods of our simulations stated in Scc. lIII Cl which does not take into account the muon-charged 
current reactions. In order to treat the muon-charged current reactions, one has to solve the neutrino-transfer for 6 
species, which requires formidable computing resources. It is also necessary to prepare the EOS table as a function 
of not only Y e but also Y^, which makes the setups further complicated. The muon-charged current reactions should 
be included in future to assess their effects. 

The second case of the treatment of muons is for so-called neutron star (NS) matter. In this case, the neutrino less 
f3 equilibrium at zero temperature is achieved and the chemical potential of the muons, fi^, is the same as that of 
electrons. Thus the relations of the chemical potentials become as follows: 

H e = fx n - (12a) 
/i M = [i e . (12b) 

Moreover, the charge neutrality condition pip is modified to 

(Y e + Y^)p B = (1 - X)YC,HPB,H + XY C ,QpB,Q, (13) 

where Y" M is a muon fraction. 

III. RESULTS 

In this section, we show the features of the EOS described in the preceding section. We choose a mixed EOS with 
the hadronic matter including pion and the quark matter of the bag constant B = 250 MeV fm~ 3 (B 1 ^ 4 = 209 MeV), 
as a reference model. This value is adopted in another recent study on the hadron-quark phase transition [34j . Here 
we discuss the properties of EOS with the applications to the maximum mass of compact stars and a stellar collapse. 
We mainly show the results of the reference model for the NS matter and the SN matter. Incidentally, results of the 
model without pions are also shown, and the dependence on the bag constant is studied for several topics. Further 
studies on stellar collapse will be reported elsewhere 0, [2(| • 

A. Assessments of EOS 

First of all, we examine the composition of our EOS for NS matter and SN matter to compare with other EOS's in 
the previous studies using similar schemes [1, H, H, H3| ■ For SN matter, we show the cases of T = 50 MeV with the 
electron type lepton fraction Yi =0.1, where Yi is defined as the sum of the electron fraction, Y e , and the electron-type 
neutrino fraction, Y Ve . Note that when the neutrinos are fully trapped, Yi is conserved for each fluid element in the 
stellar core. In FIG.[TJ we show the particle fractions, Yi = ^, where Hi represents the number density of the particle 
i, for the reference models. When pions beg in to condensate in the NS matter, the proton fraction increases and the 
fractions of electron and muon decrease 27]. At the onset of the hadron-quark phase transition, d quarks have larger 
population than those of u- and s quarks not only for the NS matter [ij but also for the SN matter [f| . In the pure 
quark phase, the number ratio of u-, d- and s quark is nearly equi-partition and the relation, Y s < Y u < Yd holds for 



2 We remark that the /3 equilibrium assumed in previous sections is for the electron- type leptons. This assumption is not inconsistent 
with the treatment of muons. 
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FIG. 1: Particle fractions for the reference model of NS matter (left) and SN matter with T — 50 MeV and Yi = 0.1 (right). 
The particle fraction, Yi, is defined as — , where m represents the number density of the particle i. 



the NS matter. For the SN matter with neutrino trapping, u quarks have larger population than those of d- and s 
quarks 0, Q. All of these features are consistent with the previous results. 

In FIG. [2J we compare the EOS for the reference model with those for pure hadronic and quark matters, where 
the bag constant B = 250 MeV fm~ 3 . This figure illustrates both for NS matter and SN matter with T = 20, 50 
and 120 MeV and Yi = 0.1. As mentioned already, the pressure of the hadron-quark mixed phase is not constant 
in an isothermal process which is well presented in our models. In addition, this trend is shown also in the 
other previous work while the bag constant is different from ours. The internal energy, which is the difference 
of the energy density with respect to pbc 2 , of the mixed phase is larger than that of the hadronic phase for high 
temperature regime (T > 50 MeV). This is not unphysical because we should compare not the internal energies but 
the free energies. In the left panels of FIG. [SJ we show the Hclmholtz free energies per baryon as a function of the 
specific volume. If the EOS is thermodynamically stable, this function is convex downward. This feature is fulfilled 
for our models. In the right panels of FIG. [3J on the other hand, we show the Gibbs free energies per baryon as a 
function of the pressure. We can recognize that the free energies of the mixed phase are always lower than those of 
the pure hadronic and quark phases. 

For the phase transition of matter with two components, it is known that the mixed phase should satisfy not only 
the condition of mechanical stability, as verified above, but also the stability with respect to diffusion [35J. In our 
model, the condition is expressed as 



< 0, (14) 



where X p is a "total" proton fraction and defined as 



y 1 ( 2n » -n d -n s \ 

X p = — \n p + j. (15) 

This feature is shown for our results in FIG. H and satisfied also in other ranges of the parameters. The reference 
model of the SN matter with Yi = 0.1 for fixed pressure and temperature are indicated in this figure. From these 
analyses, we can confirm that our EOS is thermodynamically stable. 

In FIG.0 we show the temperature dependences of the transition density and the critical baryon chemical potential 
for our EOS of the SN matter with Yi — 0.1 and 0.3. These phase diagrams are given not only for the reference model 
but also for a model without pions while the bag constant of both models is 250 MeV fm~ 3 . The baryon chemical 
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FIG. 2: EOS's for mixed matter (solid lines), pure hadronic matter (long-dashed lines) and pure quark matter (short-dashed 
lines). In all panels, the bag constant is chosen as B = 250 MeV fm -3 for the quark model. The left and right panels show the 
pressure and the internal energy, respectively, as a function of the density. In both panels, the upper left, the upper right, the 
lower left and the lower right plots correspond to the reference models of NS matter, SN matter with T = 20, 50 and 120 MeV, 
respectively. The electron type lepton fraction, Yi, is fixed to 0.1 for all models of SN matter. 



potential in these figures, is the same as the neutron chemical potential in the preceding sections, //„. We can see 
that the reference model has larger transition density and baryon chemical potential than those of the model without 
pions. This is because pions make the EOS softer; in other words, pressure gets lower for a fixed baryon number 
density. Thus the condition of equilibrium @ is satisfied for larger density. Incidentally, the phase diagram has been 
studied in Ref. [fj while its EOS for hadronic matter does not include pions. The color dielectric model is used for 
the quark matter model whereas we use the MIT bag model. In Ref. [fj, the transition density gets lower at higher 
temperature, which is the same as our result, however, its transition density is lower than ours, especially at finite 
temperature. For instance, the transition density is ~ 10 14 g cm -3 at T ~ 30 MeV. 

It is suggested that the transition line has an end point in the high temperature regime. The nature of the 
transition changes at this point, while the exact phase diagram is not well known. This is a so-called critical point. The 
temperature of the critical point, T r , is investigated experimentally by heavy- ion collisions [36|,[37| and theoretically by 
lattice QCD calculations [H, IM S3 S3 • From these studies, it may be in the range, 150 MeV <T C < 200 MeV. Above 
the critical temperature, the quark phase may be the most stable state for any densities. It is also suggested that the 
baryon chemical potential of the critical point is much smaller than the typical hadronic scale, hb 40 MeV. Although 
our model cannot describe the critical point in detail, the critical baryon chemical potential drops dramatically with 
the temperature in the regime T > 100 MeV. For much higher temperature regime, the quark matter is more stable 
than the hadronic matter even at zero density. As shown in FIG. [51 for instance, the hadron-quark phase transition 
occurs at T = 150 MeV and the quark phase is always the most stable state at T = 200 MeV for the model with 
B = 250 MeV fm~ 3 . Therefore, the plots for T > 150 MeV are not given in FIG. \5\ In this study, the bag constant 
B is independent of the temperature and the temperature at which the quark matter becomes most stable depends 
solely on B. The larger the value of B, the higher this temperature becomes. When this temperature is in the range 
150 MeV <T< 200 MeV, the bag constant is in 210 MeV fm- 3 < B < 650 MeV fin" 3 (200 MeV < B 1 / 4 < 265 MeV) 
for our models. It is noted that the bag constant of the reference model, B = 250 MeV fm -3 (B 1 ^ 4 — 209 MeV), lies 
within this range. 
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FIG. 3: EOS's for mixed matter (solid lines), pure hadronic matter (long-dashed lines) and pure quark matter (short-dashed 
lines). In all panels, the bag constant is chosen as B = 250 MeV fm -3 for the quark model. The upper left and upper right 
panels show the Helmholtz free energies per baryon (F) as a function of a specific volume of the system and the Gibbs free 
energies per baryon (G) as a function of the pressure, respectively, for the reference model of NS matter. The lower two panels 
are the same as upper two panels but for the reference model of SN matter with T — 50 MeV and Yi — 0.1. 



B. Maximum mass of the hybrid stars 

We examine the maximum mass of hybrid stars constructed by our EOS for NS matter. In the following, we denote 
the EOS's without pions and quarks (the original Shen EOS [Ty, without pions but with quarks, with pions and 
without quarks and the model with pions and quarks (the reference model) as OO, OQ, PO and PQ, respectively. 
The bag constant is B — 250 MeV fm -3 for the models with quarks. In FIG. [71 we show the mass-radius trajectories 
for our EOS's together with the recent data of compact star masses. We can see that the pion population and the 
hadron-quark transition lower the maximum mass because they soften the EOS. For instance, the maximum mass of 
model OO is 2.2 solar masses (M©) while those of models PO, OQ and PQ are 2.0Mq, 1.8M© and 1.8-Mq, respectively. 
In the right panel, we show the dependence on B. We can see that the maximum mass becomes lower when the bag 
constant becomes small. This is because the phase transition density is lower for smaller bag constants. 

The recent measurements of compact star masses increase the lower limit of the maximum mass. For instance, it 
has been established with 95% confidence that at least one of the two pulsars, Ter 5 I and J, is more massive than 



I.68M0 from analysis of the joint probabilities [42j. For the neutron star EXO 0748-676, it has been reported that the 
lower limits on the mass and radius are M > 2.10±0.28M^ and R > 13.8 ± 1.8 km with lcr error bars, while there are 
uncertainties in the analysis of the X-ray burst spectra (43J. We note that the reference model (B — 250 MeV fm -3 ) 
is consistent with these measurements whereas models with lower bag constants B < 200 MeV fm -3 in our study 
produce lower maximum masses. 



C. Application to stellar core collapse 



Recently, numerical studies on black hole formation by stellar core collapse have been done for massive stars with 
an initial mass of ~ 40M© pH I22I 0, [45| . The realistic EOS and the realistic progenitor models are used in the 
studies of spherically symmetric collapse in Refs. [2l|, [22| whereas the parametric EOS and polytropic initial model 
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FIG. 4: Profiles of neutron chemical potential, \i n , as a function of "total" proton fraction, X p , for fixed pressure P — 
220 MeV fm" 3 and temperatures T = 50 MeV (solid line) and 63 MeV (dashed line). This is a result for the reference model 
of SN matter with Y t = 0.1. 



are employed in those of axisymmetric collapse in Refs. HH, EH. Collapse of very massive ~ 300M Q first generation 
stars in the universe (so-called Population III stars) have been also studied [ly, E3, S, H§] . It is noted that these 
studies do not include the hadron-quark phase transition. In this study, we perform numerical simulations for a 
spherically symmetric stellar collapse using the EOS's examined in the preceding section. We set the bag constant 
B = 250 MeV fm -3 for all the models with quarks. A result of an evolutionary calculation for a Population III star 
with M — IOOM0 [50j is chosen as the initial model of our simulation. This model is the intermediate one between 
the models stated above. The numerical simulation of its collapse has been already done under the Shen EOS [33|. 

We follow the scheme of Ref. in our simulations. We solve t he g eneral relativistic hydrodynamics and neutrino 
transfer equations simultaneously under spherical symmetry [HJ, |52| ] Neutrino reactions are treated in detail [32] • 
Incidentally, we have assumed for the EOS that the electron-type neutrinos are in equilibrium with other particles in 
the hadron-quark mixed phase and the pure quark phase. Hence after the phase transition occurs, we do not compute 
the neutrino distribution functions and assume that they are Fermi-Dirac functions for all species conserving the 
electron-type lepton fraction, Yi, for each fluid element. Moreover, we neglect the entropy variation from the neutrino 
transport. We can justify this gap in the neutrino treatments because the density is high enough for neutrinos to be 
trapped anyway at the phase transition. It is also noted that we can compute the collapse up to the apparent horizon 
formation, which is a sufficient condition for the formation of a black hole. For the details of our simulation such as 
the initial condition, numerical methods and convergences, see Ref. [33l 

In FIG. [51 we show the time profiles of the central baryon mass density. These models have a bounce owing to 
thermal nucleons at subnuclear density (~ 1.4 x 10 14 g cm -3 ) and then recollapse to a black hole. The contribution 
of pions makes a difference at ~ 2 x 10 14 g cm~ 3 and that of quarks does so for larger density We can also see that 
the effect of quarks begins to work suddenly at the transition density whereas that of pions does gradually. This is 
because the thermal pions appear before the pion condensation. We note that the duration of neutrino emission is 
almost the same as the interval time from the bounce to the apparent horizon formation [2ll |22| . In Table U we show 
them for all models computed here. Since the EOS becomes softer owing to the contribution of pions and quarks, the 
interval time becomes shorter. Model OO (the original Shen EOS) takes 20% longer to recollapse than model PQ (the 
reference model in this study) does. Using this difference of the interval time, we may be able to probe observationally 
the EOS of hot dense matter in future. 

The total neutrino energies emitted from each model are shown in Table |TJ In this simulation, we assume that v T 
(p T ) is the same as {v ll ) 1 and the luminosities of and are almost identical. This is because they have the 
same kind of reactions and the difference of coupling constants is minor. Therefore we denote these four species as 
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FIG. 5: Phase diagrams of our EOS of the SN matter with Yi —0.1 and 0.3 for T < 150 MeV. Solid lines represent boundaries 
of hadronic matter and mixed matter and dashed lines do those of mixed matter and quark matter. For the upper plots in each 
panel, an EOS with pions is used while an EOS without pions is used for the lower plots in each panel. B = 250 MeV fm -3 is 
chosen for the bag constant for all panels. 



v x collectively in Table QJ taking the average. Comparing models OO and OQ or models PO and PQ in Table QJ we 
can see that the total energies of emitted neutrinos for the models including quarks become lower than those of the 
models without quarks. This is because the collapse is hastened by quarks, and their effects make differences only for 
the last moment as discussed above. Thus the duration of the neutrino emission becomes shorter while the neutrino 
luminosity does not change much. 

This trend can be seen in the effects of pions also. However, the neutrino luminosity for the models with pions 
differs gradually in time from that for the models without pions because the effect of pions works gradually. The 
neutrino luminosity summed over all species is given approximately by the accretion luminosity Lf, cc ~ GM V M j R v 
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FIG. 6: Same as FIG.[3]but for SN matter with T = 150 MeV and Yi = 0.1 (upper panels), and for SN matter with T = 200 MeV 
and Yi =0.1 (lower panels). 



TABLE I: Results of the numerical simulations. The definition of the models is given in the text (Sec. IHFBj) . t ICC represents 
the interval time from the bounce to the apparent horizon formation. E Vi is the total energy of emitted Vi, where E Vx = E v — 
Ev — E Vt — Ea T . Eaxi is the total energy summed over all species. 



Model 


t rcc (msec) 


(ergs) 


Ep e (ergs) 


Evx ( cr 8 s ) 


E M (ergs) 


OO 


402 


9.42 x 10 b2 


7.89 x 10 b2 


4.40 x 10" 


34.9 x 10 b2 


OQ 


352 


8.03 x 10 52 


6.58 x 10 52 


3.44 x 10 52 


28.4 x 10 52 


PO 


351 


8.05 x 10 52 


6.63 x 10 52 


3.67 x 10 52 


29.3 x 10 52 


PQ 


336 


7.60 x 10 52 


6.22 x 10 52 


3.34 x 10 52 


27.2 x 10 52 



[53], where G, R„, M and M v are the gravitational constant, the radius of the neutrino sphere, the mass accretion 
rate and the mass enclosed by i?„, respectively. We can roughly regard the neutrino sphere as the surface on which 
neutrinos are emitted, and R v is defined as, 



^mfp(r) 3' 



where R a is the stellar radius and l m f p is the mean free path of neutrino. Incidentally in our case, as is the case for 
an ordinary supernova, the proto-neutron star is formed before the black hole formation, and neutrinos are emitted 
mainly on the surface of the proto-neutron star. Thus we can regard R v and M v as the radius and mass of the 
proto-neutron star, respectively. On the other hand, the density of the proto-neutron star of the model with pions is 
greater than that of the model without pions comparing at the same time because the EOS becomes soft under the 
influence of pions as seen in FIG. [H Thus the radius of the proto-neutron star becomes smaller without changing the 
mass of it. As a result, the neutrino luminosity gets higher for the models including pions. For instance, at 300 ms 
after the bounce, the luminosity of model OO is 1.09 x 10 53 erg s _1 whereas that of model PO is 1.27 x 10 53 erg s _1 . 
While a simple comparison cannot be done because the quarks also affect the neutrino emission slightly, we can see 
that model PO has higher total energy than that of model OQ although the interval times are almost the same. In 
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FIG. 7: Mass-radius trajectories for our EOS's of NS matter. In the left panel, long-dashed, short-dashed, solid and dot-dashed 
lines represent the models OO (the original Shen EOS [l6|,[I3])' OQ, PO and PQ (the reference model), respectively, where the 
definition of the models is given in the text. In the right panel, we show the results of the models with pions, and the solid line 
is the same as that in the left panel (PO). Other lines correspond, from bottom to top, to the models with pions and quarks 
with B = 150, 200, 250, 300, 400 MeV fm -3 . In both panels, horizontal dotted lines represent the lower limit of the maximum 
mass of compact stars determined by pulsars Ter 5 I and J with 95% confidence, and the plots with lcr error bars are for the 
measurements of neutron star EXO 0748-676. 
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FIG. 8: Time profiles of the central baryon mass density. The notation of the lines is the same as in the left panel of FIG. [7] 
and the time is measured from the point at bounce. 
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conclusion, pions make the total energy of the emitted neutrinos lower by shortening the interval time and higher by 
increasing luminosity. 

The computations shown here are done for only a single model, M = IOOMq and the bag constant B = 
250 MeV fm -3 . For the models with a smaller bag constant, an interval time of the neutrino emission may be- 
come much shorter because of the smaller maximum mass of the hybrid stars. It is noted that we are preparing more 
detailed studies for various cases of stellar collapse [n| H3] . 

IV. SUMMARY 

We have constructed the EOS including the hadron-quark phase transition using Gibbs conditions for finite tem- 
perature. As for the hadronic phase, we have added pions to the EOS by Shen et al. [H, Qj]], which is based on the 
relativistic mean field theory. Our results are found to reproduce the composition in the previous studies particularly 
for the pion population in the neutrino-less f3 equilibrium state at zero temperature |27l|. As for the quark phase, we 
have adopted the MIT bag model of the deconfined 3-flavor strange quark matter [18| and the bag constant has set 
to be B = 250 MeV fm -3 for the reference model. As for the mixed phase, we have assumed that the equilibrium is 
achieved not only by the strong interactions but also by the weak interactions. Electrons and electron-type neutrinos 
are taken into account as being in weak equilibrium for the mixed phase and pure quark phase. Incidentally, we 
have treated the muons approximately. Our EOS is thermodynamically stable and exhibits qualitatively the desired 
properties of hadron-quark mixed matter, such as the temperature dependence of the transition density. 

We have studied the astrophysical implications of the calculated EOS. The most important effect of pions and 
quarks is to soften the EOS. Pions and quarks make the maximum mass of compact stars smaller. It is noted that the 
models including pions and quarks with the bag constant B < 200 MeV fm -3 give the maximum mass that is smaller 
than the recent measurements of compact stars. We have computed the collapse of a massive star with 100M Q using 
our EOS and found that the interval time from the bounce to the black hole (apparent horizon) formation becomes 
shorter by the inclusion of pions and quarks. As a result, they affect the total energy of the emitted neutrinos because 
of the shorter duration. We have also shown that the neutrino luminosity becomes higher under the effect of pions 
because they raise the density of the proto-neutron star core. We stress that we may be able to discuss the EOS of 
hot dense matter detecting the neutrinos from black hole progenitors in future. 

The current paper is the first attempt to clarify the effects due to the quark-hadron phase transition on the black 
hole formation and the neutrino emission. We are planning to study the initial mass dependence (l9| and the bag 
constant dependence [2(| of stellar core collapse systematically. We will also investigate effects of hyperons in the 
hadronic matter (2p| . There are more problems to be solved regardless of our project. For instance, other quark 
models like a color superconducting phase and effects of surface tension and screening of charged particles on the 
equilibrium condition in the mixed phase are worth investigating. 
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